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ABSTRACT 

Galaxy clusters are usually detected in blind optical surveys via suitable filtering 
methods. We present an optimal matched filter which maximizes their signal-to-noisc 
ratio by taking advantage of the knowledge we have of their intrinsic physical prop- 
erties and of the data noise properties. In this paper we restrict our application to 
galaxy magnitudes, positions and photometric rcdshifts if available, and we also ap- 
ply the filter separately to weak lensing data. The method is suitable to be naturally 
extended to a multi-band approach which could include not only additional optical 
bands but also observables with different nature such as X-rays. For each detection, 
the filter provides its significance, an estimate for the richness and for the redshift even 
if photo-z are not given. The provided analytical error estimate is tested against nu- 
merical simulations. We finally apply our method to the COSMOS field and compare 
the results with previous cluster detections obtained with different methods. Our cat- 
alogue contains 27 galaxy clusters with minimal threshold at 3-sigma level including 
both optical and weak-lensing information. 

Key words: cosmology: theory - cosmology: observations - large scale structure of 
the Universe - galaxies: clusters: general - galaxies: luminosity function, mass function 
- gravitational lensing: weak 



1 INTRODUCTION 

Clusters of galaxies lie in the densest regions of the Universe. 
In the standard ACDM model, they form from the highest 
peaks of the primordial density fluctuations. Their abun- 
dance as a function of mass and redshift is an extremely 
important tool to determine cosmological parameters and 
the history of structure formation. It is therefore fundamen- 
tal to build robust tools to detect galaxy clusters in different 
observational bands. 

In optical surveys, clusters of galaxies can be succes- 
fully identified via overdensities of galaxies with respect to 
the mean field and via the coherent distortion of background 
galaxies due to gravitational lensing. These two observables 
come together with the same data but have very different 
properties. On one hand lensing analyses have the advan- 
tage of directly mapping dark matter structures, but to be 
performed they need high-quality imaging data, and their 
results are intrinsically very noisy. On the other hand, de- 
tecting clusters through galaxy overdensities is much less 



demanding from an observational point of view, because it 
does not require any shape measurement. Moreover, if an op- 
tical cluster finding is performed on a lensing survey field, it 
can be used as a useful tool to eliminate spurious detections. 

A common approach in lensing is to convolve galaxy el- 
lipticities with a suitable filter to find the cluster signal be- 
yond the noise induced by intrinsic ellipticity of galaxies and 
cosmic shear by large-scale structures. Different filters have 
been proposed so far, including the Gaussian filter (Miyazaki 
et al. 2002), the aperture mass in the version of Schneider et 
al. (1998), its modifications by Schirmer et al. (2004), and 
the one proposed by Hennawi &: Spergel (2005). Maturi et 
al. (2005) proposed an optimal filter, which minimizes noises 
due to both the intrinsic ellipticity and the cosmic shear, se- 
lecting the scales on which the cluster signal-to-noise ratio 
is expected to be maximum. A comparison of the perfor- 
mance of the different filters used to find dark matter halos 
in cosmological simulations has been presented in Pace et 
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al. (2007). An application to GaBoDs data of the optimal 
filter is in Maturi et al. (2007). 

The idea of this work is to apply the linear filtering tech- 
nique to optical galaxy catalogues, possibly accounting for 
3D information, if available through photometric redshifts, 
to perform a combined analysis of galaxy overdensities and 
weak lensing data. 

To build an optimal linear filter for galaxy overdensities, 
we start from the pioneering approach proposed by Postman 
et al. (1996), who developed an optical filter to analyse the 
data of the Palomar Distant Cluster Survey. Many authors 
have proposed different modifications of the original Post- 
man filter: for example Kepner et al. (1999) introduced Pois- 
son (instead of the Gaussian) statistics and a rough usage of 
photometric redshifts; Dong et al. (2008) allowed the radius 
of the filter to change to adapt to the detected signal; Milk- 
eraitis et al. (2009) tuned the parameters using mock light- 
cones extracted from the Millennium Simulation (Springel 
et al. 2005). 

In this work, we construct a Postman-like filter as a lin- 
ear filter, using the estimated richness of the cluster and not 
the likelihood of the data as main object of analysis. This 
allows us to keep the linearity with respect to the galaxy 
density and to calculate an error estimate analytically. We 
also introduce the information from photometric redshifts in 
the procedure. We apply the filter to some simulated cata- 
logues to verify the consistency of the error estimates and 
to assess the reliability of the cluster redshift determination 
obtained with our algorithm, even when galaxy photometric 
redshifts are not available. The final step of this work is then 
the application of both optimal filtering algorithms (for op- 
tical and weak lensing data) to the COSMOS field (Scoville 
et al. 2007a). Previous cluster detections on the same field 
were done by Finoguenov et al. (2007) combining optical 
and X-ray information, by Olsen et al. (2007) and Scoville 
et al. (2007b) from optical data only, and by Kasliwal et al. 
(2008) and Miyazaki et al. (2007) via lensing analyses. 

The paper is organized as follows. In Section 2, we show 
the idea, the mathematical derivation and our implementa- 
tion of the filter. We also describe how to derive a redshift es- 
timate for the detected galaxy clusters, both with and with- 
out the information coming from photometric redshifts. In 
Section 3, we investigate the filter properties and test the al- 
gorithm using numerical simulations. In Section 4, we briefly 
summarize the properties of the optimal filter for weak lens- 
ing detections, highlighting the similarities with the optical 
filter. In Section 5, we apply both algorithms to galaxy cata- 
logues from the COSMOS survey and we discuss the result- 
ing cluster detections. We compare these results with previ- 
ous cluster catalogues obtained with different techniques on 
the same field. Summary, conclusions and future prospects 
of this work are provided in Section 6. 



2 OPTIMAL FILTER FOR GALAXY 
OVERDENSITIES 

To derive the optimal filter, we assume that positions, mag- 
nitudes (at least in one band), and possibly photometric 
redshifts of the galaxies in a surveyed field are available. 
Furthermore, we assume that bound structures like galaxy 
clusters present in the field can be traced by the galaxy dis- 



tributions within them and that the spatial distribution of 
background galaxies is random. Under these assumptions, 
given the characteristics of the survey, bands of observation 
and depth, it is possible to define a cluster model, describ- 
ing the expected spatial and magnitude distributions of the 
cluster galaxies and a noise model describing the background 
distribution of field galaxies. 

Thus we define a model n m (6, m) for the observed dis- 
tribution of galaxies as a function of the position and 
magnitude m, given by the sum of a field and a cluster com- 
ponent (Postman et al. 1996): 

n m (0, m) = rif(m)+An c (0, m) = n f (m)+AP(9—0 c )<f>(m) ,(1) 

where n / (m) is the magnitude distribution of field galaxies, 
C is the cluster centre, P is the projected radial profile of 
the cluster galaxies number density, <j)(m) is the cluster lumi- 
nosity function and A is a richness parameter, proportional 
to the total number of cluster galaxies. 

The detection of a cluster in the input catalogue is based 
on the comparison of the observed distribution of galaxies 
to the model n m . More specifically, for a fixed cluster centre 
C , we compute the value of A for which our model distribu- 
tion best describes the observed distribution n d (6,m). The 
likelihood of the observed data rid given the model n m is 
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assuming that the background noise is Poissonian noise and 
that the galaxy density is high enough such that the Gaus- 
sian approximation holds. 

We derive the value of the richness A that maximizes C 
by imposing 
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Only the first term in the numerator depends on the spatial 
distribution of the data with respect to 6 C . Thus, Equation 
(4) can be re- written as 



A 



= J <&(8 — C , m)rid(6, m)dQ.dm — B , 
where $ is the optical linear filter defined as 



$(0-0 c ,m) 



J nf J n f (m) 



(5) 



(6) 



and B is the contribution of background galaxies that is 
subtracted to work on a zero-mean noise field, 



dQdm . 



(7) 
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If we insert A obtained from Equation (4) back into Equation 
(2), we obtain the corresponding value of the likelihood 



C = C + 



/ — (nj, — nAdCldm 
J nj 



(8) 
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where Co is a (negative) constant that does not depend on 
the position of the cluster C . Comparing Equations (8) and 
(4), we see that the varying part of C is proportional to 
A 2 . This squared dependence of C on A reflects the fact 
that the \ 2 approach returns high likelihoods also for galaxy 
underdensities, for which A would be negative. 

This is why we consider convenient to search for galaxy 
clusters as peaks of the distribution of A as a function of 
the sky position G c , instead of likelihood maxima. In fact, 
we take full advantage of the linear response to the data, 
provided by Equation (5), and of the physical direct inter- 
pretation of A, which is a suitably normalized number of 
visible galaxies belonging to the cluster (see Section 2.1). 

The variance of A is given by <j\ = ((A — A) 2 ) , where A 
is an estimate of the true value of A, given by Equation (5). 
Since An c (m, 6) and n/(m) are random realizations of the 
cluster and field populations, the resulting a\ has two contri- 
butions: the first term is induced by the random fluctuations 
in the background, the second by the random sampling of 
the cluster model: 

-^rdQ,dm 
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-dfldm 
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-dQ.dm 



n f 



(9) 



By applying the formalism of linear matched filters, one can 
see that the filter minimizes the variance of the estimated A 
due to random fluctuations only. The second term in Equa- 
tion (9) gives the intrinsic Poissonian fluctuation of the clus- 
ter realization. 
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Figure 1. Radial profiles of the cluster model at different red- 
shifts, as labelled in the figure. The profile is flattened inside 
100 Kpc/h to avoid the central divergence, while the cut 
off at large radius corresponds to the i?200 limit. 
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Figure 2. Magnitude functions of the cluster model at different 
redshifts, as labelled in the figure. 



2.1 Modelling the galaxy clusters 

The cluster optical model is specified by the spatial and 
luminosity distribution of cluster galaxies. For simplicity, 
we assume spherical symmetry. We also assume that the 
projected radial distribution of the cluster members fol- 
lows an NFW model (Navarro et al. 1997; Bartelmann 1996; 
Meneghetti et al. 2002): 

2p s r s 



E(x) 



with 
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(10) 



: arctan 



■ arctanh 




where x = r/r s . This distribution, apart from the normaliza- 
tion p a , depends on the scale radius r s only, corresponding 
to 7?20o/c, where -R200 is the scale where the galaxy den- 
sity is 200 times the critical density, and c is the so-called 
concentration parameter. Our choice is motivated by the 
analysis made by Hansen et al. (2005) on a large sample of 



optically-selected clusters from the SDSS. They found that 
the NFW model is a good description of the cluster galaxy 
distribution up 7?200- The concentration parameters derived 
from galaxy distributions described by the NFW model de- 
pends on the richness but it is in general smaller (1 < c < 3) 
than those found in numerical simulations for dark matter 
halos of similar mass. In this context we use a scale radius 
r s — 500 kpc/h, that represents a good estimate for quite 
rich clusters (see Figs. 7 and 8 in Hansen et al. 2005). The 
resulting radial profiles computed at different redshifts are 
shown in Fig. 1, as examples. 

For the cluster luminosity distribution, we assume the 
Schechter luminosity functions (Schechter 1976), with pa- 
rameters taken from Popesso et al. (2005) , who analysed 97 
clusters observed in the SDSS. Since these are calculated 
from galaxy counts inside a radius of 1 Mpc/h from the 
cluster centre, we apply the same cut to the radial density 
profile. For a realistic concentration parameter c = 2, this 
implies that we limit our analysis within -R20O1 where the 
projected NFW is a good description of observed clusters, 
as discussed. 
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Apart from the cluster model, the other ingredient of 
our recipe is the distribution of field galaxies, that we as- 
sume to be random with a fixed mean angular number den- 
sity. The magnitude distribution of the field galaxies rif(m) 
can be estimated from the distribution of the whole galaxy 
sample, provided that the field is large enough such that the 
contribution of cluster galaxies can be considered negligible. 

2.2 Map making 

Clusters are searched in three-dimensional space, where 
their positions are given by the angular coordinates C and 
by the redshift z c . We create maps of A estimates for differ- 
ent values of the redshift, and we search for cluster detec- 
tions in these maps. Comparison of detections in different 
redshift slices is done subsequently (see Section 2.3.1). 

For every redshift, we adapt the cluster model n c in a 
proper way. First, we transform the cluster absolute lumi- 
nosity function into an apparent magnitude distribution (see 
Fig. 2), and re-scale the radial profile of the cluster model 
according to the angular-diameter distance. Then, we nor- 
malize properly the cluster model n c to set the richness pa- 
rameter A to be the total number N of visible galaxies in a 
cluster. Integrating the cluster model, we get: 



N 



= A 



A j n c (G,m)d 2 9dm 
jj P(0)d 2 9 



m)dm 



(11) 
(12) 



To set A = A, we impose the integrals of the angular dis- 
tribution and the luminosity function to be normalized to 
unity, that is 



f 11 
Jo 



p(6)d 2 e = i 



4>(m)dm 



(13) 
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where mn m is the limiting magnitude of the sample. 

At a fixed redshift z c , we build a two-dimensional map 
evaluating A over a grid of angular positions 9 C . Since we 
deal with discrete quantities, the integral in Equation (4) 
must be approximated as a sum over the galaxy positions. 
Only the first term in the numerator needs to be evaluated 
for every grid position 6 C , while the second one is equal to 
unity because of the normalization. The denominator for a 
given redshift is a constant for all sky positions. We denote 
it by C. Thus the discretized version of Equation (4) reduces 
to 

A* = A{9i,9 s ) = ^± , (15) 



where 



Mi - 



n c (r k ,m k ) 
n f (m k ) 
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n f (m k ) 



(16) 



In the formula above, r k is the angular distance of the galaxy 
from the centre 



r k = \0 k - {9 h 9j)\ 



(17) 



and the index k runs over all N galaxies of the sample. Clus- 
ters candidates are identified with the peaks of the resulting 
A-map. To select only significant peaks, we calculate their 
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Figure 3. Analytic response of the algorithm (Eq. 18) for two 
clusters at redshifts 0.2 and 0.4. When the cluster is seen at the 
correct redshift, the value of the response is unity by definition 
(black points). 



signal-to-noise ratios using the variance given in Equation 
(9). 



2.3 Redshift determination 

2.3.1 Without photometric redshifts 

To estimate the cluster redshift, we cannot just maximize 
A with respect to z c . Indeed, this would introduce a bias as 
discussed below. The analytic response of the algorithm as 
a function of the search redshift z c for a cluster at a redshift 
z c and with a galaxy distribution n c is 



A(z c ) = 



— n c d£ldm 

Thf 



(18) 



<>f 



-dQdm 



This corresponds to Equation (4) when applied to a known 
distribution rid = rif + n c . By construction, A = 1 when z c = 
Zc, but there is no reason why this should be a maximum for 
the function A(z c ). This is shown in Fig. 3, where we plot 
the analytic response for two clusters located at z = 0.2 and 
z = 0.4. For both of them the maximum of A as a function 
of z c is found at a higher redshift than the cluster's. 

For a correct redshift estimate, the likelihood £ (Equa- 
tion 8) in the peaks of the 2D distribution of A must be used. 
By definition, the value of £ is maximal when the cluster is 
filtered by a model at the correct redshift (See Fig. 4). 



2.3.2 With photometric redshifts 

If the photometric redshift is available for each galaxy, we 
can use this further information in our analysis. Assuming 
that a redshift probability distribution p k (z) is derived for 
the fc-th galaxy of the sample, p k (z c ) can be used as a weight 
factor for that galaxy when estimating the field number den- 
sity, rif, the cluster richness, A, and the likelihood, £, at a 
redshift z c .. 

On the other hand, this can be interpreted as an ex- 
tension of the two-dimensional linear matched filter to a 
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dfldmdz 
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(24) 



In the process of estimating the redshift from multi- 
band photometry, a spectral type is also assigned to each 
galaxy. This can be used to k-correct the values of the mag- 
nitude and to obtain rest-frame values, which are directly 
comparable to the model. The filter is applied to the data 
as usual, creating maps for different redshift slices and se- 
lecting peaks, as described in Section (2.2). The most likely 
redshift estimate for a detection is then choosen. 



Figure 4. Analytic value of C for two clusters at redshift 0.2 and 
0.4, as a function of the search redshift z c - 



three-dimensional one, where the redshift probability func- 
tion represents the cluster profile along the line of sight. 
In fact, in an ideal case where all galaxies have the same 
redshift probability distribution, which for simplicity we as- 
sume to be a Gaussian distribution centered on their most 
probable location, z e3 t,h, with a constant rms, a z , 



Pk{z) = —j=, — exp 
v2na z 



{Zest,k — z) 2 
2<7? 



(19) 



the three-dimensional model with a cluster at redshift z c is 

n m (0,m,z) — nf(m,z)+An c (9,m,z) 

= n f (m,z)+AP(0-0 c )<t>z(m)p(z), (20) 

where 

(Zc - zf ' 
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If we then further assume that the background distribution 
is only slowly changing with redshift, so that the background 
n/(m, z) for z c — a z < z < z c + a z is constant, the richness 
estimate at redshift z c is 



A : 
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n/(m, z c ) 



(rid{6, m, z) -— nj(m, z c ))dQdmdz 
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This is equivalent to Equation (4), but with the additional 
dimension. In practical applications, Equation (22) has to 
be discretized in analogy with Equation (15). Thus, 



Ai 



= E 



P(r k )4>(m k )p k (z c 
nf(m k , z c ) 



(23) 



In this case the numerator is evaluated weighting every 
galaxy with p k (z c ). The constant denominator C, that does 
not depend on the observed galaxy spatial distribution, can 
be evaluated using a model for cr z , possibly depending on 
the magnitude and the redshift estimate of the galaxies. The 
same model is also necessary for the evaluation of the vari- 
ance that, in analogy with Equation 9, is 



2.4 Comparison with other optical cluster finding 
methods 

The optimal matched filter we used for cluster detections 
(both optical and weak lensing) relies on the general knowl- 
edge we have of galaxy clusters, as incorporated in the clus- 
ter model, and selects all physical properties which enables 
us to distinguish them from the field. It is in fact a general 
approach which can incorporate at once different cluster as- 
pects in contrast to other methods which aim at a single 
signature such as galaxy overdensities, cluster galaxy red- 
sequence, cD galaxies, etc. In other words, instead of en- 
forcing some criteria which is later applied to data, we first 
define the general properties of galaxy clusters and then let 
the filter find what are the distinctive features for their de- 
tection given the actual data noise. 

In the following sections we discuss the main techniques 
applied in literature and how they compare with the one 
presented in this work. 



2.4-1 Galaxy overdensities 

Looking for galaxy overdensities in astronomical images 
by eye was the first method used to detect galaxy clus- 
ters. In the last decades different automated methods were 
proposed, such as friend-of-friend algorithms (Li & Yee 
2008), Voronoi tesselation (Ramella et al. 2001), filtering by 
wavelets (Finoguenov et al. 2007), adapting kernels (Mazure 
et al. 2007) or matched filters (e.g. Postman et al. 1996; Kep- 
ner et al. 1999; Gilbank et al. 2004; Dong et al. 2008; Menan- 
teau et al. 2009; Milkeraitis et al. 2009). Together with the 
additional features discussed in the following sections, our 
algorithm incorporates the mentioned matched filters. In 
fact it includes the angular clustering information on the sky 
according to the P term of Equation (23) or, if photometric 
redshifts are available, the full three-dimensional clustering 
information thanks to the p term of Equation (23). In more 
details, the main differences with respect to the matched 
filters mentioned in this section are: 

(i) our filter can take advantage of the information from 
photometric redshifts in a very flexible way, accounting for 
their uncertainty. Instead of simply slicing the data-set and 
neglecting all galaxies outside a fixed redshift range (e.g. 
Milkeraitis et al. 2009), we consider all galaxies belonging 
to the data set by weighing them according to their own 
redshift probability distribution. 
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(ii) It uses distinct luminosity functions for galaxy clus- 
ters and field galaxies to define the cluster model and the 
noise contribution. Although this seems an obvious choice, 
previous filtering algorithms use the luminosity functions 
taken from the whole population of galaxies, including the 
field, to define the cluster model. 

(iii) It is linear with respect to the data, because it makes 
use of the richness A instead of the likelihood C as the main 
object of analysis. In fact, even if they are based on the same 
statistal distribution, A is linear, in contrast with C which 
is quadratic, and therefore is capable to distinguish positive 
overdensities from negative underdensities to which C would 
assign high probabilities. 

(iv) It gives as a natural output an estimate of the num- 
ber of observed galaxies in a cluster, which can be easily 
corrected for redshift dimming to get a meaningful physical 
quantity (see Section 5.1). 

2.4.2 Brighter cluster galaxy 

The algorithm used by Koester et al. (2007), based mainly 
on red-sequence information, includes a cut-off term for 
the brighter cluster galaxy luminosity, such that detections 
without a very brilliant central galaxy are discarded. In our 
approach, instead of looking for a specific cluster member, 
we let the filter select in an optimal way the whole luminosity 
function 4> expected for a cluster at a given redshift in con- 
trast with the field luminosity function n/, as in Equation 
(23). As the most massive and luminous galaxies are found 
preferably at the center of the clusters, at low magnitudes 
the ratio between the cluster model luminosity function and 
the observed field distribution will increase. The algorithm 
will assign large weights to very luminous galaxies and thus 
they will be a strong indication of the presence of clusters. 
The advantage of our approach is that the filter is defined 
according to well defined statistical quantities and not to 
a free parameter. In addition we allow the algorihm to be 
sensitive to a broad range of systems, from rich clusters to 
groups which might otherwise fall below a given threshold. 

2.4-3 Cluster galaxies red sequence 

Red-sequence cluster finders rely on the fact that the cen- 
ter of galaxy clusters is typically populated by elliptical 
galaxies whose color is on average redder than the one of 
the field galaxies at the same redshift. Usually a cut in a 
colour-magnitude diagram is performed to select these red 
galaxies overdensities that identify the galaxy cluster candi- 
dates. This technique, presented in Gladders & Yee (2000), 
has been succesfully used in different galaxy surveys (e.g. 
Gladders & Yee 2005; Koester et al. 2007; Lu et al. 2009; 
Thanjavur et al. 2009). We thus mention and explain how 
this technique can be implemented in the framework of our 
method with a simple generalization, even if we did not ex- 
plicitly made use of color information in this work. In a 
following paper we will fully discuss a multi-band approach 
which deserves an extended discussion on its own. 

Our method belongs to the matched filters developed for 
CMB observations (e.g. Sch&fer et al. 2001; Pace et al. 2008) 
which use at once all available bands to obtain a unique re- 
sponse. With this respect, the algorithm presented here is a 



special case where a single band is considered, but could 
be naturally extended to include the full available color 
information. In this way, we would consider not only the 
two bands that are usually used in classical red-sequence 
searches, but all bands and color properties of the clusters 
and the field, adding new possible information. Again, this 
would be done not by defining a cut-off in the color space 
but by defining the proper weight to each galaxy performing 
an optimized analysis. 

In any case, already at this stage without multi-band 
information, we partially account for the massive galaxies, 
as those targeted more precisely by the red sequence, thanks 
to the use of the cluster luminosity function (see Section 
2.4.2). 



3 INVESTIGATING THE FILTER 
PROPERTIES 

The final step of this work will be the application of the 
algorithm to real data, namely the i' band galaxy catalogue 
of the COSMOS field (Ilbert et al. 2009), limited at i' = 25. 
Before doing this, we have to evaluate the capability of our 
algorithm to detect clusters and to correctly measure their 
richness and redshift. In these tests we use the field distribu- 
tions measured from the COSMOS data to define the filter. 
Thus, we use the specific implementation of the algorithm 
that will be later adopted for the final data analysis. 

3.1 Signal-to-noise estimates 

We now evaluate the expected signal-to-noise ratios for clus- 
ters with different redshift and richness, as described by 
our cluster model, once the actual properties of the galax- 
ies in the COSMOS catalogue are assumed. The expected 
signal-to-noise ratios when the photometric redshifts are ne- 
glected are shown in the upper panel of Fig. 5. In this case, 
the expected signal and noise are given by Equations (4) 
and (9), respectively, and the signal-to-noise ratio peaks at 
z ~ 0.45. To explain this behaviour, we plot separately in 
Fig. 6 the expected signal of a cluster with 100 visible galax- 
ies at z = 0.2 and the two noise terms as a function of red- 
shift. The signal, i.e. the number of visible galaxies of the 
modeled cluster, decreases for increasing z as expected (solid 
line). The trend for the noise is different: on one hand, the 
term related to the background galaxies only (dot-dashed 
line) grows monotonically as a function of redshift up to 
z ~ 0.5 after which it basically stays constant; on the other 
hand, the intrinsic fluctuations of the cluster signal decreases 
monotonically with z. Asa result the latter term is dominant 
at low redshifts (up to 2 < 0.4). 

This noise behaviour derives from the fact that 
for low-redshift clusters, the filter, being proportional to 
n c (m)/n/(m), selects the bright end of the cluster luminos- 
ity function, that has almost no background. Therefore the 
background noise decreases with redshift while the intrinsic 
fluctuations of the cluster signal increase because the A es- 
timate depends in practice on the very few bright galaxies 
only. 

In order to consider the case in which photometric red- 
shift information is available, we use the values for the red- 
shift (zest) and its 68% confidence levels (z m in, z max ) re- 
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ported in the COSMOS catalogue to associate to each galaxy 
a probability distribution given by 



p(z) 



where 



and 



2Tva„ 



exp 



(z — Z e st)' 

2a 2 



Zest Zmin if 2 <^ Zest 
Z-max ~ %est l£ Z Z es f 



= (z„ 



i in )/2 



(25) 



(26) 



(27) 



When estimating the normalization constants, i.e. the 
denominator of Equation (22) and the variance (Equation 
24), which do not depend on the observed spatial galaxy dis- 
tribution, we use the mean redshift error for different classes 
of objects, estimated from the comparison with a spectro- 
scopic subsample (Ilbert et al. 2009). For our dataset, these 
errors are: 



= < 



0.007 for 0.2 < z < 1.5, 

0.011 for 0.2 < z < 1.5, 

0.053 for 0.2 < z < 1.5, 

0.06 for z > 1.5 



i' < 22.5 
22.5 < i' < 24 
24 < i' < 25 



.(28) 



As we can see in the lower panel of Fig. 5, using the 
photo-z information, the picture becomes more confused, 
as the sensitivity depends on the field population density 
nf(m,z) at the search redshift (see Equation 24). The only 
clear trend is the decrement of the sensitivity at high red- 
shift due to the dimming of the cluster signal. In this case 
we do not have any worsening of the performance at low 
redshift because the weight given by the filter through the 
n c (m)/rif(m, z) ratio selects almost the same part of the 
luminosity function at each redshift. 

3.2 Simulations with mock catalogues 

We then test the application of our algorithm to mock galaxy 
catalogues, built from the COSMOS data. To create realistic 
catalogues for the field galaxies, we first randomize the posi- 
tions of the galaxies to cancel any structures contained in the 
field. We assume that the cluster galaxies do not affect the 
photometric properties of the overall sample because their 
number is negligible with respect to that of the field galaxies. 
On the top of this random background, we add some mock 
clusters. The cluster galaxies are distributed according to 
the model for the spatial and luminosity distribution of the 
cluster members (see Section 2.1). The photometric redshifts 
are assigned following a Gaussian redshift probability distri- 
bution with errors given by Equation (28). For illustrative 
purposes, we first generate a square field with side 1 degree, 
containing 25 galaxy clusters of different richness and red- 
shift placed on a regular grid. In Fig. 7 we show the results 
of the analysis of this field without the usage of photometric 
redshifts, at a search redshift z c = 0.4. All clusters leave an 
imprint on the A-map, whose strength depends on the rich- 
ness and the redshift, with clusters at z — z c being brighter 
than the others with the same A. 

Other simulations are done to verify the linear response 
of the filter with respect to the cluster richness. In this case, 
we assume for simplicity all clusters at z = 0.4 and com- 
pare the input number of galaxies to the estimated one. The 
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Figure 5. Expected signal-to-noise values for the detection of a 
cluster, without (top panel) and with (bottom panel) photometric 
redshifts, for the COSMOS i' band catalogue. The black contours 
refer to different integer values for S/N, starting from unity at the 
bottom (see the color bar on the right). The value on the y-axis 
is the number of galaxies of the cluster that are detected when it 
is located at z = 0.2. 



results are shown in Fig. 8 where we also verify the good 
agreement of the analytic estimate of the variance given in 
Equation (24) with the one resulting from our Monte Carlo 
simulations. 

We also use simulations to test the capability of the al- 
gorithm to estimate correctly the redshift of a cluster. The 
results are shown in Fig. 9. The mock clusters were com- 
posed requiring that they have 200 galaxies under the mag- 
nitude limit at z — 0.2. At higher redshifts, fewer galaxies 
are seen and analysed, explaining the increasing variance of 
the measurement. Note that even with single-band observa- 
tions a redshift estimate for the clusters is possible albeit 
with larger errors with respect to the case in which the pho- 
tometric redshift information is included. A similar result 
was obtained by Dietrich et al. (2007). 
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Figure 6. Signal and noise for a cluster with 100 visible galaxies 
at z = 0.2 as a function of redshift. The signal (in red) corre- 
sponds to the number of galaxies below i' = 25, calculated from 
the model luminosity function. The dot-dashed line represents 
the noise due to the background galaxies, while the dashed line 
represents the noise due to the fluctuations in the cluster galaxy 
distribution. The solid blue line represents the total noise. Scales 
for noise and signal are shown on the left and right, respectively. 
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Figure 7. Distribution of the richness A at z c = 0.4. Colorbar 
for A is reported on the right. The grid identifies the input posi- 
tions of the mock clusters. The labels along the axes indicate the 
richness and redshift of the cluster at that position. 



4 OPTIMAL FILTER FOR WEAK LENSING 

In the next Section we will apply to the COSMOS field both 
optimal filters for galaxy cluster detection from optical and 
weak lensing data. Here we introduce weak gravitational 
lensing and then the filtering method we will adopt for lens- 
ing data. 

The observable quantity in weak lensing analyses is 
the reduced shear g, a measure of the distortion of back- 
ground galaxies by intervening structures (i.e. galaxy clus- 
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Figure 8. Calibration of the values for the richness A values ob- 
tained by the algorithm with and without redshift information. 
Clusters are located at z = 0.4. The shaded area represents the 
analytic estimate of the uncertainty, when the photometric red- 
shift information is used. 
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Figure 9. Redshift determination obtained by the algorithm, 
using the likelihood C instead of the richness A, with and without 
redshift information. When photometric redshift are used, error 
bars become negligible. 



ters). From the 2D projection of the cluster mass distribu- 
tion, £(#), one can calculate the lensing potential 



4G DiD, 



d 2 6>'E(0')hi|0-0'| , 



(29) 



c 2 D ls 

where Di, D s and Di a are the angular diameter distances 
between lens and observer, source and observer, and lens 
and source, respectively. The complex reduced shear is then 
defined as 



9(0) 



7 (0) 



<0) 



(30) 



where the shear 7 and the convergence k are second-order 
derivatives of the lensing potential, namely 



K = -{(pll 



71 = o (0U ~ 022) , 



72 = <f>12 



(31) 

(32) 
(33) 
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using the notation 



dOidOj 



(34) 



The linear matched filter for weak lensing cluster detec- 
tions has been presented in Maturi et al. (2005). We refer 
the reader to that paper and to Pace et al. (2007) for a de- 
tailed description of its specific derivation, the comparison 
with other filtering techniques and tests of its performances. 
Here we summarise the main properties of the filter and how 
it can be applied to real data. 

The weak lensing filter for cluster detection is expressed 
in Fourier space for convenience. Its shape is given by 



*(fe) 



1 



(27T) 



P N (k) 



P N (k) 



(35) 



where f is the Fourier transform of the weak lensing signal 
expected from a lensing cluster, i.e. the reduced shear, and 
Pjv represents the power spectrum of the noise, due to the 
intrinsic ellipticity of the sources, their finite number on the 
sky and the additional shear induced by large-scale struc- 
tures. The weak lensing signal model for clusters is computed 
assuming a spherically symmetric NFW density profile for 
the dark matter distribution (Equation 10). Again, as in the 
case of galaxy overdensities, the estimate of the lensing sig- 
nal is obtained by convolving the shear data, D(0), with the 
filter *, 



A{0) = j D{0')^(0 -0')d 2 6' . 
The variance of the estimate is 



2 



1 

(2^ 



/ l*(fc)| : 



P N {k)d 2 k 



(36) 



(37) 



The filter in Equation (35) was constructed requiring that 
the estimate of the amplitude of the signal A given by Equa- 
tion (36) is unbiased and has minimum variance. The optical 
filter is the solution of the same problem for a galaxy distri- 
bution in magnitude and position. Note that if we neglect 
the noise contribution given by large-scale structures, i.e. we 
assume the noise to be white, this lensing filter and the one 
for the galaxy overdensities are fully equivalent, except of 
course for the specific nature of the two. 

As in the case of the optical filter, the application of the 
weak lensing filter to the data is done evaluating A(0) over 
a grid of angular positions. Equation (36) is then discretized 
to be applied to real data. It becomes 



1 N 

A(9) = -J2^(\^-0\) , 



(38) 



where n g is the number density of galaxies used for weak 
lensing measurement and e t k denotes the tangential compo- 
nent of the fe-th galaxy ellipticity with respect to the angular 
position. Thus, the filtering is performed in the real domain, 
which requires to back-Fourier transform the filter function 
defined in Equation (35). For every estimate of A, its er- 
ror a\ is also computed, according to Equation (37). Peaks 
with S/N > 3 are selected from the map and considered as 
detections. 



5 COSMOS FIELD 

We now apply the optical filter and the weak lensing filter to 
the COSMOS field, a 2 square degree equatorial field that 
has been observed with the Hubble Space Telescope and 
with many other instruments, covering wavelengths from X- 
ray to radio (Scoville et al. 2007a). 



5.1 Optical detections 

For the optical analysis, we use the public galaxy catalogue 
with photometric redshifts (Ilbert et al. 2009), considering 
the angular position, the photometric redshift, and the i'- 
band magnitude measured with Subaru (Taniguchi et al. 
2008). The photometric redshifts presented in the catalogue 
were obtained from the galaxy fluxes in 30 different bands, 
ranging from the ultraviolet to the near infrared. To each 
galaxy a SED has been assigned during the determination of 
the photometric redshift. We use this additional information 
to calculate the k-correction. We search for clusters between 
z — OA and z = 0.8, with steps of Az — 0.02. The upper 
limit of z = 0.8 is motivated by two factors: first, the model 
that we assume for the cluster luminosity function is based 
on low-redshift objects, and cannot be extrapolated to too 
high redshifts; second, if the k-correction becomes too strong 
for a given type of galaxies, our sample can no longer be 
assumed complete up to a given rest-frame magnitude. We 
scan an area smaller than the full field of view, such that, for 
each redshift, the filter (cut at r — 1 Mpc//i) is completely 
inside the survey field. In this way we avoid any border effect, 
at the price of loosing some possible detections of structures 
that have their centre inside the COSMOS area, but which 
are not completely contained in the field of view. 

In Fig. 10, as examples, we show the results for the slices 
at z = 0.48, 0.50, 0.52 where the detections with S/N > 3 are 
marked with circles. Note that some structures are visible 
in more than one slice. In that case we assign as redshift 
of the detection the one where the likelihood is maximal, as 
discussed is Section 2.3. In this way we obtain 140 significant 
detections with a redshift estimate. 

Since the absolute magnitude limit changes with red- 
shift, the richness of clusters with the same galactic popu- 
lation is redshift dependent. In particular the richness de- 
creases as the redshift increases. To correct for this effect and 
to obtain a quantity which depends on the cluster galaxy 
content only, we make use of the luminosity function (j>( m ) 
of a cluster at an arbitrary redshift z c = 0.2 to compute the 
normalization factor 



R(z) = 



m lim( z ) 



4>(m)dm . 



(39) 



This quantity represents the part of <f>(m) visible at a given 
redshift z. Dividing the measured richness A for R(z), where 
z is the estimated redshift of the cluster, we obtain a value 
that is directly proportional to the physical galactic content 
of a cluster. More precisely this 'corrected richness' corre- 
sponds to the number of cluster galaxies that would be vis- 
ible if the cluster was located at z = 0.2. The 'corrected 
richness' of our detections is shown in Fig. 11, together with 
the selection threshold S/N = 3 for optical detections al- 
ready derived and shown in the lower panel of Fig. 5. Notice 
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Figure 11. Redshift and (corrected) richness of the detected clus- 
ters. The black line represents the detection limit, that depends 
on model and field distributions at different redshifts. Clusters 
found in the optical analysis only are marked in red, those with 
a counterpart in the lensing analysis are marked in blue. 



T 


| 1 


I 


i 


i 1 





i r 

X 




X 








+ 




- x 


+ 

: * 


+ + 

+ x + 


+ 

- +* 


- 


■ 5- 




' + > 

\ + 


< *: 

+ 




+ 

++ +■ 




it +x 




+ I 















that the hills and wells of the curve are due to the change 
of the field population as a function of redshift. 



Figure 10. Maps of the estimates of the richness A on the COS- 
MOS field, at redshifts z = 0.48, 0.5, 0.52, from top to bottom. 
Red circles indicate significant peaks, while red crosses represent 
a cluster detection at that redshift. 



5.2 Lensing detections 

The weak lensing cluster detections in COSMOS are per- 
formed by applying Equation (38) as discussed by Maturi et 
al. (2005) on the shear catalogue obtained with the Suprime- 
Cam mounted on SUBARU (Miyazaki et al. 2007). The 
galaxy number density in the catalogue is 42 arcmin -2 . Fol- 
lowing Miyazaki et al. (2007), we assume that the mean red- 
shift of the sources is 0.8. The weak lensing cluster model 
used in the filter construction is the signal expected from a 
cluster of mass 5 x 1O 14 M0 at z — 0.5. 

The total number of significant (e.g. with S/N > 3) de- 
tections is 82. The expected number of spurious detections 
at this S/N level is around 40%, as computed by Pace et 
al. (2007). The derived weak lensing detections are related 
to the optical overdensities using a correlation length of 500 
kpc//i. This quite severe constraint is justified by the ab- 
sence of any redshift information in the lensing data. The 
27 detections satisfying this correlation criterion are listed 
in Table 1, while optical detections without a lensing coun- 
terpart are listed in the Appendix. 

In Fig. 12 we show the redshift distribution of our com- 
bined detections. From the source redshift distribution it 
is possible to evaluate the redshift sensitivity function for 
weak lensing, that is a measurement of how much the shear 
of observed galaxies is sensitive to structures at different 
redshifts. In Fig. 12 we also show the optical sensitivity 
function, computed as the S/N of a cluster with 100 visible 
galaxies at z = 0.2. The ratio of lensing-confirmed clusters 
over the total number of optical detections in our analysis 
does not show any clear trend with redshift, but this could 
be due to the poor statistics. 



5.3 Comparison with literature 

We can compare our results to published catalogues of clus- 
ters previously extracted from the COSMOS field. Previous 
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Figure 12. Top panel: redshift distribution of all optical detec- 
tions (dashed line) and of the lensing-confirmed clusters (solid 
line). Bottom panel: lensing sensitivity function and optical sen- 
sitivity function, defined as the S/N of a cluster with 100 visible 
galaxies at z = 0.2, in arbitrary units. 

analyses were done by Finoguenov et al. (2007), combining 
optical and X-ray information and by Olsen et al. (2007), 
who found galaxy overdensities, by applying a Postman- 
like filter to a previous version of the galaxy catalogue with 
no redshift information. Scoville et al. (2007b) analysed the 
COSMOS field with a wavelet method, and presented a cata- 
logue of large-scale structures. We note that these structures 
are in general much larger than our targets, and that their 
size in some cases exceeds 10 Mpc in comoving coordinates. 
The few previous cluster detections in this field via lensing 
were reported by Kasliwal et al. (2008), while comparing the 
performance of ground-based and space-based telescopes for 
weak lensing, and by Miyazaki et al. (2007). 

Correlating our detections with results from literature, 
we use different criteria depending on the characteristics of 
the catalogues. For sources presented in Finoguenov et al. 
(2007) and Olsen et al. (2007), we match detections using a 
physical correlation length of 1 Mpc//i, i.e. the cut-off ap- 
plied to our radial filter. Instead, for the larger structures 
presented in Scoville et al. (2007b), we use a correlation 
length equal to the FWHM of their density peaks. In all 
these cases, the correlation along the line-of-sight is per- 



Table 1. Catalogue of the clusters detected with both optical 
and lensing filters. The angular position is the one of the detection 
from optical data. In the last column we indicate whether the clus- 
ter has been found previously with other techniques: 'X' stands 
for Finoguenov et al. (2007), 'S' for Scoville et al. (2007b), 'O' for 
Olsen et al. (2007), 'L' for Kasliwal et al. (2008) or Miyazaki et 
al. (2007). 



formed considering differences in redshift smaller than 0.05. 
This can be quite severe for Olsen et al. (2007), as their 
redshift determination was based on single band observa- 
tions without any redshift information. For the lensing de- 
tections presented by Kasliwal et al. (2008) and Miyazaki et 
al. (2007), we use the same criterion we adopted for the in- 
ternal matching with our lensing detections, i.e. a correlation 
length of 500 kpc//i. Regarding Finoguenov et al. (2007), we 
must stress that the comparison can only be done with their 
catalogue of X-ray confirmed clusters, obtained correlating 
420 optical overdensities with 150 X-ray diffuse sources. This 
justifies some discrepancies between the two analyses. 

The existence of previous detections of the objects 
found in our optical and lensing analysis is reported in the 
last columns of Tables 1 and Al. 

We note here the ability of our optical filter to distin- 
guish different galaxy clusters that are aligned along the 
line of sight, despite their signals are degenerate in weak 
lensing analyses. In fact, the two detections with identifi- 
cation numbers 4 and 7 correspond to SLJ1000. 7+0137, a 
weak lensing and X-ray source studied by Hamana et al. 
(2009) with spectroscopical observations. They found that 
it is actually made by two, or possibly three, overdensities 
of galaxies at different redshifts, almost at the same angular 
position. Our method is able to disentagle the two structures 
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Figure 13. S/N distribution of all the significant optical detec- 
tions (solid red), and of the detections with a counterpart in other 
analyses (dotted red). In blue we show the same distributions, 
but considering only combined (optical+lensing) detections. The 
black crosses indicate the percentage of our optical detections 
having at least one counterpart (scale on the right). 



Figure 14. Mass distribution of all the X-ray detected clusters 
(in black), of the clusters we detect in the optical analysis (in 
red), and of those we detect in both optical and lensing analyses 
(in blue). The black crosses indicates the percentage of clusters 
detected in our optical analysis. 



at z = 0.22 and z = 0.34 using information from photomet- 
ric redshifts only. A similar case is that of SLJ1001. 2+0135, 
described by Hamana et al. (2009) as the superposition of 
two overdensities, located at z — 0.22 and z = 0.37. Our de- 
tection 8 corresponds to the latter, while at lower redshifts 
that angular position is outside the scanned area, because 
the filter would not be completely contained in the field of 
view (see Section 5.1). However, if we force our algorithm to 
analyse that region, we find an overdensity at z — 0.24 with 
coordinates (R.A. = 150.329, Dec. = 1.609), confirming the 
alignment of two different structures. 

5.4 Characterizing our sample 

In Fig. 13 we plot the number of detections obtained with our 
analysis (both optical and optical plus lensing, dashed and 
solid lines respectively), as a function of their optical S/N 
ratio together with the corresponding amount of detections 
already reported in literature. As expected, the percentage 
of 'matched' detections increases as a function of S/N. The 
only detection with S/N > 5 without a counterpart is ac- 
tually near to a large structure observed by Scoville et al. 
(2007b) with a redshift mismatch of Az — 0.06, only slightly 
larger than the limit we adopted. 

To investigate the selection function of our method 
with respect to the mass, we plot in Fig. 14 the number 
of detections with respect to the X-ray masses estimated 
by Finoguenov et al. (2007). We include in this analysis 
only those X-ray clusters which lie inside the volume of our 
search, avoiding borders of the field and too high redshifts, 
as already discussed in Section 5.1. We see that the rate of 
X-ray clusters that are optically confirmed by our analysis 
increases with their mass, and it is above 50% for clusters 
of mass above 1.5 x 10 13 M Q . 

Finally we compare in Fig. 15 the corrected richness, 
corresponding to a measurement of the mass in galaxies, 
with the X-ray mass estimates made by Finoguenov et al. 




1 10 

X mass (10 13 M Q ) 



Figure 15. Correlation between the X-ray mass derived by 
Finoguenov et al. (2007) and our corrected richness. In red we 
show all the optical detections, in blue the lensing-confirmcd clus- 
ters. 



(2007). The proportionality between the two 'mass esti- 
mates' is clear although the scatter is very large. Note that 
the corrected richness is only a relative measure of the clus- 
ter stellar mass, in the sense that we do not quantify the 
physical mass scale. Note also that our filter has a fixed 
physical size, while the X-ray analysis was done inside an 
estimated r$oo for each cluster, partially explaining the dif- 
ferences. The noise in the richness estimate for low-mass 
clusters is likely to be due to the field galaxies we observe 
inside our fixed spatial filter. This is expected given the error 
bars that show our analytic estimate. 
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6 SUMMARY AND CONCLUSIONS 

In this paper we have presented an optimal linear filtering 
technique for the detection of galaxy clusters from optical 
and weak lensing data. The filter was first presented by Ma- 
turi et al. (2005) for weak lensing analysis and we now use 
the same formalism to extend it to optical data. The filter 
relies on a physical model for clusters and thus it accounts 
for all the known properties of galaxy clusters which can 
be inferred from photometric observations. Different infor- 
mation can be included and in this work we decided to re- 
strict our analysis to the following observed properties of the 
galaxies: positions, magnitudes in one band, and photomet- 
ric redshifts if present. The algorithm starts from a galaxy 
catalogue with the mentioned physical quantities and cre- 
ates a map of values for the richness parameter A at differ- 
ent redshifts. The peaks of these maps represent the possible 
locations of galaxy clusters. The redshift of the clusters is 
determined by maximizing the likelihood of the data as a 
function of the redshift of the cluster model. 

Our algorithm offers a number of improvements over ex- 
isting filtering methods. Photometric redshifts are included 
in a flexible way, that adapts itself to the precision of the 
measurement of each galaxy, avoiding any sharp cut of the 
catalogue in redshift slices. It uses results from observed 
galaxy clusters to build a suitable model for cluster galaxy 
distributions. The richness calculated by the algorithm at 
the position of the cluster is a measurement of the number 
of galaxies in the detected cluster, that can be corrected for 
redshift dimming. Moreover, from the expected spatial pro- 
file and luminosity distribution of the clusters, and from the 
observed field population, it is possible to evaluate the noise 
of the detection and its significance. We tested the algorithm 
with numerical simulations, probing that it is able to obtain 
an unbiased estimate of the richness of the cluster, with an 
uncertainty predicted analytically from the model, and also 
its redshift. 

We applied both the weak lensing and the optical fil- 
ters to real data, obtaining a catalogue of candidate galaxy 
clusters for the COSMOS field. We presented a catalogue 
of 27 lensing-confirmed clusters, 11 of which do not have 
any previous detection in the literature. For the subsam- 
ple of objects that have been analysed through their X-ray 
emission, we found that we are able to detect more than 
50% of the clusters with an X-ray mass over 1.5 x 10 13 Mq 
and we found a good correlation between our galaxy rich- 
ness parameter and the cluster mass determined from X-ray 
temperature. 



APPENDIX A: OPTICAL-ONLY DETECTIONS 

In Table Al we list the optical detections without lensing 
counterpart. 
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ID z R.A. Decl. S/N richness previous 

detections 



28 


0.12 


"i a a o c i 

149.862 


1.765 


4.279 


69.35 


X 


29 


0.12 


1 rr\ a 11 

15U.432 


2.630 


3.925 


59.13 




30 


0.16 


150.387 


2.069 


3.246 


44.64 




31 


0.22 


150.102 


2.358 


5.042 


144.29 


adU 


32 


0.22 


150.395 


2.460 


4.397 


111.20 


U 


33 


0.22 


150. 620 


2.707 


4.070 


96.12 


n 
O 


34 


0.26 


1 ACX O A 1 

149.842 


2.682 


3.436 


56.19 


O 


35 


0.26 


150. osy 


2.474 


3.640 


62.42 




36 


0.26 


150.268 


2.672 


3.305 


52.37 


c 
O 


37 


0.28 


149.594 


2.150 


3.978 


75.93 




38 


0.28 


i a a n a n 

149.969 


2.450 


3.190 


50.52 


b 


39 


0.28 


1 rr\ a A A 

15U.U44 


2.225 


4.159 


82.51 


bO 


40 


0.28 


150.194 


1.747 


3.280 


53.16 


b 


A 1 

41 


0.30 


150.186 


2.794 


3.517 


74.18 




42 


0.30 


150.572 


1.942 


3.337 


67.24 


X 


1 Q 

43 


0.32 


i rn /ice 


2.049 


5.332 


182.39 


Q 
O 


44 


0.32 


150.652 


2.313 


3.531 


83.32 


c 
b 


45 


0.34 


149.596 


2.820 


4.773 


146.29 


Xb 


46 


0.34 


1 A Ci TC"~7 

149.767 


2.329 


3.336 


75.73 


o 
O 


47 


0.34 


i p» a c\ r 7Ci 

150. 070 


2.378 


3.372 


77.19 


X 


48 


0.34 


150.185 


1.765 


4.804 


148.11 


Ab 


49 


0.34 


"1 P" Ci 1 'T O 

150.373 


2.444 


3.691 


90.87 


xs 


50 


0.34 


150.536 


2.730 


4.186 


114.54 


n 



51 


0.36 


149.787 


2.167 


3.900 


98.05 


n 
O 


•52 


0.36 


150.119 


2.689 


5.227 


168.45 


X 


•53 


0.36 


150.523 


2.570 


4.385 


121.46 




54 


0.38 


149.675 


2.413 


3.133 


66.70 


c 
O 


55 


0.38 


"1 A Ci f ~7C* r T 

149.767 


1.625 


3.112 


65.94 




56 


0.38 


i ac\ on 

149.821 


2.275 


4.637 


135.17 


bU 


•57 


0.38 


"1 A Ci 1 \ f ' / ' 

149.966 


1.678 


4.250 


115.08 


X 


58 


0.38 


150.234 


2.474 


3.609 


85.62 




•59 


0.38 


150.387 


2.413 


4.162 


110.80 


X 


60 


0.38 


150.663 


2.559 


3.077 


64.64 




61 


0.40 


1 A Ci 

149.578 


2.533 


3.343 


61.52 




62 


0.40 


1 A Ct C\ r 7-\ 

149.971 


2.748 


3.401 


63.40 




63 


0.40 


150.356 


2.651 


4.816 


119.47 




64 


0.42 


1 A Ci /~"~7 A 

149.674 


2.731 


3.336 


72.20 




65 


0.44 


149.961 


2.210 


3.315 


66.64 


xo 


66 


0.44 


1 rr\ AC\ A 

15U.494 


2.070 


5.034 


142.59 


Ab(J 


67 


0.44 


"i p" a r' act 

150.690 


2.021 


3.695 


80.85 




68 


0.46 


1 I Ci C r"~7 

149.667 


1.625 


4.658 


118.42 


b 


69 


0.46 


150.354 


2.741 


3.178 


58.86 




70 


0.46 


150.653 


2.271 


4.031 


90.50 




71 


0.46 


150.721 


2.244 


3.160 


58.26 




72 


0.48 


149.561 


2.519 


3.407 


74.82 




73 


0.48 


149.707 


2.506 


4.519 


125.21 


c 



74 


0.50 


149.525 


2.441 


3.519 


78.18 


b 


75 


0.50 


149.695 


2.657 


3.148 


64.27 




76 


0.50 


150.113 


2.559 


4.855 


140.51 


Ab 


77 


0.50 


1 PT Ci 11 O 

150.328 


2.735 


3.562 


79.87 


b 


78 


0.52 


1 A n C A o 

149.542 


1.719 


3.733 


73.37 




79 


0.52 


149.714 


2.267 


3.885 


78.72 


o 


80 


0.52 


149.823 


1.821 


3.785 


75.19 




81 


0.52 


150.059 


1.629 


3.646 


70.42 




82 


0.52 


150.141 


1.597 


4.453 


100.43 




83 


0.52 


150.295 


1.680 


5.809 


163.77 




84 


0.54 


150.135 


1.853 


4.970 


135.34 




85 


0.54 


150.216 


1.821 


4.928 


133.21 


X 


86 


0.54 


150.260 


1.765 


3.528 


72.84 




87 


0.54 


150.360 


1.627 


3.883 


86.33 




88 


0.54 


150.467 


2.066 


5.088 


141.32 


X 


89 


0.54 


150.573 


2.166 


3.591 


75.14 
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ID 2 R.A. Decl. S/N richness previous 

detections 



90 


0.56 


149.525 


1.760 


4.173 


123.00 




91 


0.58 


149.578 


1.683 


3.028 


78.39 


S 


92 


0.60 


149.617 


1.740 


5.134 


184.96 


S 


93 


0.60 


150.063 


2.798 


3.874 


110.18 




94 


0.60 


150.253 


2.346 


3.334 


84.34 




95 


0.60 


150.301 


2.804 


3.445 


89.34 




96 


0.60 


150.491 


2.745 


3.851 


109.03 


X 


97 


0.60 


150.574 


2.608 


3.827 


107.78 




98 


0.60 


150.616 


2.780 


3.527 


93.14 




99 


0.60 


150.717 


2.537 


3.799 


106.39 




100 


0.62 


150.052 


2.321 


3.082 


77.73 




101 


0.62 


150.520 


2.473 


3.620 


104.23 




102 


0.62 


150.707 


2.759 


3.306 


88.25 




103 


0.62 


150.725 


2.625 


3.492 


97.54 




104 


0.64 


150.633 


2.715 


3.738 


100.25 




105 


0.64 


150.737 


2.825 


3.644 


95.81 


X 


106 


0.66 


149.802 


1.804 


3.204 


83.66 


s 


107 


0.66 


150.185 


2.158 


3.438 


94.33 




108 


0.66 


150.196 


2.238 


3.613 


102.74 




109 


0.66 


150.687 


1.667 


3.469 


95.77 




110 


0.68 


149.673 


2.277 


3.670 


92.24 




111 


0.68 


149.925 


2.597 


3.092 


70.21 


s 


112 


0.68 


149.948 


2.097 


4.162 


113.53 




113 


0.68 


150.010 


2.119 


3.810 


98.07 




114 


0.68 


150.060 


2.608 


6.305 


234.42 


so 


115 


0.68 


150.173 


2.518 


5.706 


195.95 


s 


116 


0.68 


150.257 


1.968 


3.895 


101.70 




117 


0.70 


149.647 


2.828 


3.514 


93.08 




118 


0.70 


149.964 


2.673 


4.606 


146.79 


s 


119 


0.70 


149.986 


2.578 


5.200 


181.54 


XSO 


120 


0.70 


150.003 


2.451 


5.684 


212.81 


s 


121 


0.70 


150.152 


2.601 


3.509 


92.83 


s 


122 


0.72 


149.899 


2.394 


3.789 


118.18 


s 


123 


0.72 


150.086 


2.460 


4.479 


157.29 


s 


124 


0.72 


150.108 


2.565 


3.475 


102.36 


s 


125 


0.72 


150.207 


2.361 


3.056 


83.12 


s 


126 


0.72 


150.593 


2.129 


3.011 


81.17 




127 


0.72 


150.736 


2.416 


3.731 


115.19 


s 


128 


0.74 


149.523 


2.656 


3.344 


101.46 




129 


0.74 


149.866 


2.492 


4.466 


170.82 


s 


130 


0.74 


149.986 


2.563 


3.305 


99.40 


so 


131 


0.74 


150.024 


2.688 


3.388 


103.83 


s 


132 


0.74 


150.051 


2.302 


3.821 


128.61 


s 


133 


0.74 


150.116 


2.705 


4.030 


141.57 


so 


134 


0.74 


150.171 


1.703 


3.449 


107.13 




135 


0.74 


150.296 


2.378 


3.196 


93.72 


s 


136 


0.74 


150.547 


2.803 


3.212 


94.55 


s 


137 


0.76 


149.559 


1.647 


3.054 


88.26 




138 


0.76 


149.732 


2.761 


3.248 


98.66 




139 


0.76 


149.839 


1.684 


3.468 


111.29 




140 


0.76 


150.114 


2.255 


3.080 


89.62 


s 


141 


0.80 


149.706 


2.265 


3.251 


91.47 




142 


0.80 


150.040 


2.652 


3.906 


126.31 




143 


0.80 


150.368 


2.005 


3.546 


106.42 


X 


144 


0.80 


150.506 


2.222 


3.229 


90.37 




145 


0.80 


150.538 


2.148 


3.097 


84.12 


o 


146 


0.80 


150.580 


2.652 


4.594 


169.49 


s 


147 


0.80 


150.702 


2.769 


3.850 


123.09 





Table Al. As Table 1, but for clusters detected with the optical 
filter without lensing counterpart. 
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